Electron Self-injection in Multidimensional Relativistic Plasma Wakefields 
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We present an analytical model for electron self-injection in nonlinear, multidimensional plasma 
wave excited by short laser pulse in the bubble regime or by short electron beam in the blowout 
regime. In this regimes, which are typical for electron acceleration in the last experiments, the laser 
radiation pressure or the electron beam charge pushes out background plasma electrons forming a 
plasma cavity - bubble - with a huge ion charge. The plasma electrons can be trapped in the bubble 
and accelerated by the plasma wakefields up to very high energies. The model predicts the condition 
for electron trapping and the trapping cross section in terms of the bubble radius and the bubble 
velocity. The obtained results are in a good agreement with results of 3D PIC simulations. 

PACS numbers: 52.38.Kd,52.65.Rr,52.27.Ny 



Plasma-based charged particle acceleration is now a 
fast developing area of science. It has attracted much 
attention due to the recent experimental breakthrough 
in generation of quasimonoenergetic, dense and short 
bunches of relativistic electrons with up to GeV energies 
in the laser-driven acceleration [l], 0,13, 0] and due to 
the energy doubling of a 42 GeV electron beam [3] . The 
high-gradient acceleration is provided by very strong elec- 
tromagnetic fields in a plasma wake excited by short laser 
pulse or electron beam. The plasma electrons are radi- 
ally expelled by ponderomotive force of the laser pulse 
or by Lorentz force from the space charge of the electron 
beam that leads to formation of the plasma cavity - bub- 
ble with unshielded ions inside. This is the bubble @ 
or blowout 0] regime of laser-plasma and beam-plasma 
interaction. It is important that a small part of the back- 
ground plasma electrons can be trapped in the bubble 
providing the electron self-injection. The bubble velocity 
is close to the speed of light and the trapped relativistic 
electrons can be continuously accelerated in the plasma 
fields up to a very high energy. 

The electron self-injection is the key phenomenon of 
the last laser-plasma acceleration experiments [H, [3, 13, 0] . 
It is a crucial factor for the quality of the accelerated 
electron beam. Many applications ranging from x-ray 
free electron lasers to electron-positron colliders require 
high quality electron beams with a very low emittance 
and energy spread. The one-dimensional model based on 
the Hamiltonian formalism has been developed to study 
electron trapping in relativistic plasma wave 

UM- Sa- 
cral optical and plasma techniques have been proposed 

to enhance the electron self-injection in plasma wakefield 
like the collision of two counter-pr opag ating laser pulses 
T5L 11 1, density transitions 13, 13, 14 1 etc. In the ultra- 
high intensity regime the radial force from the driver and 
from the plasma fields strongly affect electron motion and 



the electron dynamics become complex and multidimen- 
sional. The trajectory of the trapped electrons starts at 
the bubble front, bends around the cavity and becomes 
caught at the trailing edge of the bubble. Evidently, such 
trajectories cannot be described in the framework of a 
one-dimensional model. 

Despite the great interest in plasma-based electron 
acceleration there is little theory of self-injection. In 
Ref. [H M, [l7| a multidimensional model of a rela- 
tivistic plasma wake including beam loading is proposed, 
however the electron self-injection was not considered. 
The electron self-injection has been mainly studied nu- 
merically 18, [13, Eo]. It was found in Ref. [13] that a 
sufficiently large bubble, R > 4, can trap plasma elec- 
trons. Here R is the bubble radius normalized to c/lu p , 
ujp = y/ 'Aim^e 2 / m is the (non-relativistic) plasma fre- 
quency for the background electron density no, e and 
m, are the charge and mass of the electron, respectively. 
A different condition for electron self-injection has been 
proposed in Ref. [13] : R > 7o, where 70 is the bubble 
gamma-factor. Unfortunately, principal limitations of a 
numerical approach make the accuracy and the validity 
range of the obtained conditions unclear and do not pro- 
vide an insight in the self-injection physics. However, the 
numerical approach has its natural limitations in accu- 
racy and validity conditions of the obtained results and 
does not provide insight in the self-injection physics. 

In this paper, we present an analytical model for self- 
injection and check it by direct measurements of model 
parameters in 3D particle-in-cell (PIC) simulations. For 
simplicity we assume that the bubble has a spherical 
shape with large radius R 3> 1 that is typical for relativis- 
tic laser pulses. The numerical simulations [23] and theo- 
retical analysis [13] demonstrate that the bubble shape is 
close to the spherical one for relativistically intense laser 
pulse and when the number of the trapped electrons is 
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not too large. It can be shown FjJ, l20| that the space 
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time distribution of the electromagnetic field inside the 
spherical bubble expressed in cylindrical coordinates is 
E x = (l + V)£/4, E ± = -Bg = r ± /4, where £ = x - Vt, 
and V ~ 1 is the bubble velocity. We use dimensionless 
units, normalizing the time to w" 1 , the lengths to c/l 
the velocity to c, the electromagnetic fields to mcw p /\e\ 
and the electron density, n, to no. The distribution of 
the electromagnetic field inside the plasma cavity can 
be also expressed through vector and scalar potential 
[20l | A x = r 2 /8, = with gauge A x = —ip, where 
r 2 = £ 2 + y 2 + z 2 is the distance to the bubble center. 

We assume that the electron trajectory is plane (z = 0) 
because of the axial symmetry of the wakefield and the 
driver. The driver field is neglected because it is weak 
at the region where trapping occurs (far behind the laser 
pulse or electron beam). The field of the trapped parti- 
cles is also neglected since we assume that the number 
of the trapped particle is not too large. The electron 
dynamics in the bubble is governed by the Hamiltonian 
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where II is the canonical momentum of the electron and 
£, y are coordinates canonically conjugated to II. As the 
Hamiltonian is not the function of time then it is integral 
of motion H = const. The Hamiltonian equations are: 



dp x 
dt 
d£ 
dt 
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■Py) 1 ^ 2 is the relativistic 7 
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factor 

A(£, y) is the electron 



where 7 = (l+p% 
of the electron, p = V7 = II 
momentum. 

We assume that the electron is located initially at the 
bubble border p = 0, y = R and £ = at t = 0. It 
follows from the initial condition that 



H = 1 -V Px -{l + V)v{t,y) = 1 + (1 + V) 



R 2 



(3) 



where ip(r = 00) = cp(r = R) = — R 2 /8. We solve the 
electron motion equations numerically. To be more re- 
alistic and in accordance with PIC simulations we in- 
clude an electron sheath around the plasma cavity, which 
screens the bubble ion field in the surrounding plasma. 
We model the electromagnetic fields inside the bubble 
as follows E x = /(r)£/2, E y = —H z = f(r)y/A, where 
f(r) = [tanh(i?/d - r/d) - l]/2, d is the width of the 
electron sheath. The similar distribution of the wake 
field with d < 1 is observed in the 3D PIC simulation 
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20 ]. The typical trajectory of a trapped electron is 
shown in Fig. Q] for R = 7, d = 0.3 and 7o = 4. The 
sheath does not affect the electron trajectory in the bub- 
ble when d < 2. We can conclude from the numerical 
solutions that the electron very soon becomes ultrarela- 
tivistic so that p x 3> p y 2> 1. The electron undergoes 




FIG. 1: Trajectory of the trapped (solid line) and untrapped 
electron (dashed line) calculated by numerical solution of 
equation of electron motion and the bubble border (dashed 
circle). The coordinates are given in c/uj p . 



betatron oscillations about the £-axis and slowly moves 
along this axis (\d£/dt\ — \p x /j — V\ <C 1). The condi- 
tion d£/dt > can be considered as a necessary (but not 
sufficient) one for electron trapping. The electron leaves 
the bubble and is lost when r > R. The most critical 
instants of time when the electron can leave the bubble 
is t = t m > when dy/dt = p y = and the electron 
excursion from the £— axis reaches its maximum. At this 
moment we can write 7 — Vp x ~ p x j (279) +1/ (2p x ), 
where p x » 1 is assumed, 70 = (1 — V 2 )^ 1 ^ 2 S> 1 is the 
gamma-factor of the bubble. The square of the distance 
of the electron from the bubble center (£ = 0, y = 0) at 
t — t 7r)i is 



R z 



7o 2 



Px (im ) 



(4) 



where we use Eq. (J2]). 

We change the variables p y = R 2 P y , p x = R 2 P X , £ = 
XR, y = YR and t = Rs. As a result Eqs. (P) and © 
takes a form 



dP x 
ds 

ds 

dX 
ds 




IP 2 + P 2 



+ O(i?- 4 ,7 - 2 ), (5) 



(6) 
(7) 



where P x = P y = and dP y /ds = —1/4, dP x /ds = at 
s = 0. In the zeroth order in 7~ 2 and R~ 4 Eqs. ©-© 
do not depend on any parameters and can be numerically 
solved to find s m and t m = Rs m + O (7q~ 2 , i?~ 2 ): Si — 
3.2, s 2 ^ 10.7, s 3 ~ 23.0, .... Using Eq. © the condition 
that electron leaves the bubble (r m > i?) can be written 
as follows 



2> R2Pxi 2 Sm) +0(% 4 ,R- 2 ) 
7o 



(9) 



The electron leaves the bubble if there exists at least 
one value of s m when the condition ([9]) is satisfied. We 
can take only P x (si) ~ 1.1 since P x (si) < P x (s 2 ) < 
Px (S3) < .... As a result we come to the condition for 
electron capture in the bubble 



To < J_ 

R ~ y/2' 



(10) 



which is close to the condition obtained numerically in 
Ref. [ij. 

Now we estimate the effect of the impact parameter on 
the trapping condition. If p < R then electron first moves 
through the decelerating bubble field where £ > 0. As a 
result the electron gains the negative momentum, A, at 
£ = in contrast to the electron with p = R, which starts 
motion from position £ = with P = 0. The deceleration 
leads to reduction of the longitudinal momentum at s = 
si for electron with p < R as compared to the electron 
with p = R. We estimate this reduction, P x (si,p = R) — 
P x (si, p < R), as A, that is close to the value calculated 
by the numerical integration of Eqs. |1}, ©• To calculate 
A we can use Eq. §5§ and assume l/fi< E-p< fi. As 
p y ~ at f = 0, y ~ p then H ~ (1 + i? 4 A 2 ) 1/2 + i? 2 A + 
p 2/4 = l + j R2/4 5 so that A = jy(^ + 2/i? 2 )/(2^ + 2/i? 2 ) ~ 
i//2, where 1/ = (1 - p 2 /R 2 )/A. Therefore, P x (si,p < 
R) ~ P x (si, p = i?) — v/2 and the condition for electron 
trapping (fTTj|) can be rewritten as follows 



Po < P < R, 



(11) 

(12) 



Therefore the trapping cross-section a near the trapping 
threshold 70 — 2 _1 / 2 i? takes the form 



ttR 2 



2p_ 

R 2 



dp= 1 



R 2 



'R 2 



(13) 



It follows from Eq. (|T3|) that the trapping cross-section 
decreases as jo/R increases and the trapping stops at the 
threshold 70 ~ 2~ 1 / 2 i?. 

The developed model abstracts from some important 
effects observed in numerical simulations: the bubble 
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FIG. 2: (a) The distribution of the electron density in the 
plane x — z is calculated by 3D PIC simulation, (b) The 
number of the trapped electron per laser period as function of 
xq. (c) The current coordinate, x, of the trapped electrons as 
function of xo and function approximating x(xo) (red dashed 
line), (d) The gamma- factor of the bubble as function of t 
estimated as 70 — (—2dx/dxo)~ 1 ^ 2 (red line) and the trapping 
condition 70 = 2 _1//2 i? (blue line). The coordinates and time 
are given in the laser wavelength and laser periods, the density 
is given in the critical plasma density, respectively. 
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shape deformation during laser pulse propagation and 
the bubble field enhancement at the bubble back due to 
electron sheet crossing [H It follows from 3D PIC 



simulations that the bubble back typically moves slowly 
than the bubble front nana. The bubble back gamma- 
factor should be used in Eq. (flU|) as 70 since the electrons 
are trapped at the bubble back. Another effect is that the 
bubble field can be stronger at the bubble back than it 
follows from linear approximation for the bubble field in 
our model. To estimate the effect of field enhancement on 
the trapping we introduce the field enhancement factor 
k so that A x = —tp = kr 2 /8. It follows from the scala- 
bility analysis of Eqs. ©-JE]) that the coefficient at RHS 
of Eq. (flQ|) becomes \fkj2. The considered effects can 
significantly affect electron self-injection. For parameters 
of simulations performed in Ref. [1 91 ] the minimal radius 
of the bubble with self-injection is reduced by a factor of 
4 when the bubble field enhancement (~ 4, see Fig. 3 in 
Ref. [3]) and the bubble back gamma- factor (~ 9, see 
Fig. 1 in Ref. [3]) are taken into account in Eq. (Ti"0|) . The 
estimated radius is about 1.5 times more than observed 
in the simulations. Further studies are needed to include 
the above mentioned effects in our model accurately. 

We carry out a numerical simulation of laser-plasma 
interaction by PIC VLPL3D code. The laser pulse 
is circularly polarized and has the envelope a = 
ao exp (— r 2 /r|) cos (Trf/IY), where Tl = 30/s is the 
pulse duration, — 9fim is the focused spot size and 
ao = eA/mc 2 — 1.5 is the normalized vector poten- 
tial, which corresponds to the laser intensity / = 3 x 
10 18 W/cm 2 , A = 0.82um is the laser wavelength. The 
plasma density is no = 1.16 x 10 19 cm -3 so that the pa- 
rameters are close to that used in experiments Q. The 
typical distribution of the electron density in the bubble 
regime is shown in Fig. [2]^a) when laser pulse has passed 
the distance 2000A in plasma. Because of the large num- 
ber of trapped electrons the field enhancement effect is 
suppressed and the bubble shape deviates from the ideal 
spherical one. We modify the code to track individual 
particles over the full simulation. To separate trapped 
electrons, we introduce an energy threshold: we assume 
that an electron was trapped if its maximum energy ex- 
ceeds 75 MeV. 

The distribution of the number of the electrons trapped 
in the bubble as a function of their initial position x$ is 
shown in Fig. E^b) for three instants of time. Fig. [2](b) 
demonstrates that the capture process is non-uniform 
and the number of trapped electrons changes along the 
laser path. The trapped electrons start motion from po- 
sition x = xq at t = io that corresponds to the electron 
coordinates y « R, £ = and the bubble center position 
x = xq according to our model. If the bubble propa- 
gates with velocity V(t) then xq — J Q ° V(t')dt' . For the 

particle position we have x = xq + f t V x (to, t')dt', where 
v x is longitudinal component of the electron velocity and 



v x (to,to) — 0. Differentiating this equation with respect 
to t we get dx/dxo = 1 + [1/F(t )] J* a [dv x (t ,t')/dt ]dt', 
where we use dxo/dto = V(to)- In our model V x (to,t) is 
a function of t — to and for a highly relativistic electron 
7 70 1, v x ~ 1 we can write a dx/dxo ~ — l/(27g). 
This derivation does not take into account the bubble 
shape deformation. However, it follows from our simula- 
tions that 70 calculated from dx/dxo — —1/(2jq) largely 
determines the gamma-factor of the bubble back because 
the trapped electrons with 1 < 7 < 70 spent most of 
time in the bubble back. It is seen from Fig. [2jb)-(d) 
that the self- inject ion is strongly correlated with this 
7— factor. The number of trapped particles per laser 
period peaks in the regions llOOA/c < t < 1400A/c 
and t > 1700A/c, where 70 is minimal. Vice versa, for 
1400A/c < t < 1700A/c when 7o > 2~ 1 / 2 .R the number 
of the trapped particles becomes negligible that agrees 
with predictions of the model. 

In Conclusion we present the model for electron self- 
injection in the relativistic plasma wakefield generated 
by a short laser pulse or by an electron beam. The den- 
sity threshold for electron self-injection is predicted by 
our model that is different from the self-injection condi- 
tion proposed in Ref. [IH . The evidences of such thresh- 
old have been observed in experiments llj. The model 



can also explain electron self-injection scheme based on 
downward plasma density transition 14| . Such transition 
leads to the bubble elongation because of plasma density 
decreasing Further studies are neededand sharply reduces 
the gamma-factor of the bubble back thereby strongly 
enhancing self-injection. Recently, a strong impact of 
the wake phase velocity on electron trapping was also 
observed experimentally [22} that is consistent with our 
model. However, the accurate measurements of the bub- 
ble velocity and the structure of the electromagnetic field 
in the bubble are needed for careful verification of the ob- 
tained results. Much effort is now mounted in theory and 
in experiments to control the emittance and the energy 
spread of the accelerated electron bunch. It follows from 
the obtained result that the self-injection strongly de- 
pends on the wake phase velocity that can be controlled, 
for example, by plasma density profiling [23l |. Therefore, 
the obtained results can be used to optimize the plasma- 
based acceleration. 
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